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Abstract: We examine the effects of invisible particle emission in conjunction with 
QCD initial state radiation (ISR) on quantities designed to probe the mass scale 
of new physics at hadron colliders, which involve longitudinal as well as transverse 
final-state momenta. This is an extension of our previous treatment [1] of the effects 
of ISR on global inclusive variables. We present resummed results on the visible 
invariant mass distribution and compare them to parton-level Monte Carlo results 
for top quark and gluino pair-production at the LHC There is good agreement as 
long as the visible pseudorapidity interval is large enough (?7max ^ 3). The effect of 
invisible particle emission is small in the case of top pair production but substantial 
for gluino pair production. This is due mainly to the larger mass of the intermediate 
particles in gluino decay (squarks rather than W-bosons). We also show Monte Carlo 
modelling of the effects of hadronization and the underlying event. The effect of the 
underlying event is large but may be approximately universal. 
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1. Introduction 

Amongst the many observables that could be of use in the search for new physics 
at hadron colliders, those that do not depend on detailed hypotheses about the 
structure of the final state may be the most suitable for an initial general survey of 
the high-energy frontier. Quantities of this type have been investigated and named 
global inclusive observables in ref. [2]. Such quantities can provide information on 
the energy scales of new hard processes, for example the production of exotic heavy 
particles, without imposing prejudices about the nature of the new dynamics. 

One may distinguish between global inclusive observables that depend only on 
transverse momenta, such as the visible and missing transverse energies, and observ- 
ables such as the visible energy and invariant mass, which depend also on longitudinal 
momenta. The former are not affected by the unknown motion of the hard process in 
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the collider frame, in the approximation that the process is initiated by partons mov- 
ing coUinearly with the beams. However, the colour charges of those partons and the 
high scale of the hard process necessarily imply that there is significant initial-state 
QCD radiation (ISR), which contributes to both types of global inclusive observables. 

In ref. [2] various global inclusive observables were investigated, including those 
that make use of longitudinal as well as transverse momentum components. The 
quantities studied included the visible energy E, i.e. the sum of energies of particles 
registered by the detector, and the visible invariant mass M, 



M= JE^-P^-^^, (1.1) 



where is the visible longitudinal momentum. Here ^t, the missing transverse 
energy, is in practice defined as the negative of the total visible transverse momentum. 



so that ( |1 . 1| ) does indeed define a Lorentz-invariant quantity. In addition, in [2] a 



new variable was introduced, defined as 



5mr„(Minv) = ^JM^+ + ^ML+ ^1 , (1.2) 

where the parameter Mjnv is a variable estimating the sum of masses of all invisible 
particles in the event: 

^inv 

Mi„, = ^m, . (1.3) 

i=l 

It was argued that the peak in the distribution of Sj^^^ is a good indicator of the mass 
scale of new physics processes involving heavy particle production. 

In ref. [1] we examined the effects of QCD initial state radiation (ISR) on 
these observables, first in an approximate fixed-order treatment, taking into account 
collinear-enhanced terms, and then in an all-orders resummation of such terms. We 
quantified the way the distributions of quantities that involve longitudinal momenta 
depend on the scale of the underlying hard subprocess and on the properties of the 
detector, in particular the maximum visible pseudorapidity r^max- 

In the treatment in ref. [1] we supposed that all the final-state particles from 
the hard subprocess are detected. We considered in detail the case of top quark pair 
production when both the t and i decay hadronically. In the present paper we study 
the effect of invisible particle emission, which allows us to extend the validity of our 
treatment to many cases of Standard Model and beyond Standard Model processes 
of interest at the LHC. We examine semi-leptonic/hadronic and fully semi-leptonic 
decays of the top quark (see Fig. |l|). We also consider gluino pair-production, where 
the final decay products contain heavy invisible particles, the lightest supersymmetric 
particles (LSPs), which remain undetected. The emission of the LSPs in the decay 
of the gluino can cause a sizeable effect on global inclusive observables. We examine 
the interplay of this effect with ISR effects. 
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Figure 1: Feynman diagram for qq' + ISR g ^ tt ^ jjllS/T -The solid lines represent 
visible jets or leptons, the dotted lines invisible neutrinos. 



The Monte Carlo results presented in refs. [1,2] show that the second term on the 
right-hand side of eq. ( |1.2|) is not strongly affected by ISR. The first term is intended 
to add extra longitudinal information about the hard subprocess, allowing a more 
reliable determination of its mass scale. The extra information enters through the 
visible mass M. Therefore, as in ref. [1], we concentrate on this quantity. 

The paper is organised in the following way: In section ^ we review the ISR 
resummation. Then in section ^ we develop our treatment of ISR effects to in- 
clude one or two invisible decays originating from the hard process. In section ^ we 
present comparisons of resummed distributions to distributions produced using the 
general-purpose Monte Carlo event generator Herwig++ [3] for tt and gg production. 
Moreover, we examine the effect of hadronization of ISR and the effect of the un- 
derlying event on the parton-level distributions. Our conclusions are summarized in 
section ^. Appendix ^ contains details of the numerical Mellin inversion method 
used to obtain the evolution kernels and appendix ^ gives the parton-level cross 
sections for the processes under investigation. 



2. ISR effects without invisible particle emission 

We first review briefly the main results of ref. [1]. Consider the emission of ISR 
partons from incoming partons a and b, as in Fig. ^ If all the products of the 
hard subprocess and all the ISR emitted at angles greater than 9c are detected, the 
resummed differential cross section may be written as: 

dcT f 



2.1 



3 



or equivalent ly 



dcT f 

^ dMHY ^ / ^^l'^^2-f^«'«(^l)/a(^l)Qc)-f^fe'b(22)/6(S2,<5c)<5-a'fe'(^1^2M^) , (2.2) 

where M and Y are the invariant mass and rapidity of the detected system (the 
'visible' mass and rapidity) and 

M y _ M _y . . 

Xi = —;=e , Xo = —;=e , (2.3) 

\/S being the overall cm. energy- squared. The hard subprocess cross section aa'b' 
is evaluated at cm. energy squared = X1X2S = ziZ2M'^, where Zi = Xi/xi. The 
parton distribution functions (PDFs) fa^b are evaluated at the lower scale Qc ~ OcQ} 
The kernel functions Ka'a and K^b describe the evolution of the PDFs from scale 
Qc to Q. They satisfy an evolution equation like that of the parton distributions 
themselves: 

Q^Kb,b{z) = / '^4Pb>a{z')K,b{z/z') , (2.4) 

OQ TC J Z' 

with the initial condition that Kab{z) = SabS{l — z) at Q = Qc- We describe in 
appendix ^ the method that we used to compute the kernel functions. 

The main conclusions from the study in ref. [1] are that, in the absence of invisi- 
ble particles, the above analytical results are in good agreement with those of Monte 



Carlo simulations, and that the distribution of the new variable ( |1 . 21) is indeed de- 
termined primarily by that of the visible mass M. 



3. ISR effects including invisible particle emission 

Suppose now that an invisible 4-momentum is emitted from the hard subprocess. 
If we define the total lab-frame 4-momentum of the incoming partons a and b as 

= + X2), 0, 0, (xi - X2)] , (3.1) 

then the visible 4-momentum will be ~ Pinv definition, the visible mass is 
then given by: 

= (P - p_)2 = P^P^ + - 2pl^P^ . (3.2) 

Equation ( p.2| ) demonstrates the interplay between two effects: on one hand ISR 
increases the 'true' scale of the hard process Q, to the 'apparent' scale M by con- 
taminating the detector with extra particles, and on the other hand the invisible 

^We find that results are somewhat sensitive to the constant of proportionality here. We actually 
use Qc — Q exp(— r/„iax) where ?7niax = — lntan(0c/2) is the maximum visible pseudorapidity. 
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particle emission decreases M by the loss of particles. In the case of gluino pair- 
production both effects are equally important, as we will show. 

Substituting from eq. ( |3.1| ) in eq. ( p.2| ) and defining = p?^^ ± p?^^ we obtain 

= X1X2S + m\^^ - ^fS\xvp'^^^ + ^2P+ J . (3-3) 

where rriinv represents the total invariant mass of the invisibles, ral^^ = PinvPinv,fi- 

The momenta are defined in the lab frame, relative to which the cm. frame 
of the hard subprocess is boosted by an amount defined by the momentum fractions 
xi and X2 of the partons entering the subprocess. This implies that the transform 
as: 

+ _ + - _ 1X2 _ 

Pinv A/ ™ ^inv ' Pinv A; ™ ^inv ' ) 
\ X2 \ Xi 

where = q^^v ^ifnv^ defined in terms of the invisible momentum, q^^^, in the cm. 
frame of the hard subprocess. Substituting the expressions of eq. ( p.4|) in eq. (|3.3| ) 
we find an expression for the visible invariant mass: 

= ^Inv + X1X2S [1 - Zif^^ - Z2f^nv] , (3-5) 

where we have defined = q^^/Q and used = X1X2Z1Z2S. We may now solve 
eq. (|375| ) for Q"^ to obtain Q"^ in terms of M^: 

Q = -rz — 7^^: — T~ ■ (3.6) 

The above expression for the hard subprocess scale now becomes the argument of 
the parton-level cross section, aa'h' in eq. (|2.2|): 



^77^7^= / dz^dZ2K,:a{Zl)fa{Xl^Qc)Ky,{z2)h{x2,Qc)^a'y { '''^^^l "^T^ ) 
dM^dY J \l-ZifZ,^- Z2finv J 

(3.7) 

The functions f^^^ which are related to the invisible particle four-momenta, remain 
to be determined. The visible system rapidity, Y , is also modified by the presence 
of invisible particles as: 

Y = -\og( MlziiL^^ , (3.8) 



2 Va;2(l - ^2/i„J 



and therefore eqs. (|2.3|) for 0:1,2 become 



Xl 



Y 



(M2-mL)(l-^2/,nJ 
'S'(l — ^l/mt) ^2/.m^)(l — ^l/in*;) 

^ ; - ^'inv)(^ - Zlftnv) ^-Y (^3 g^. 
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The kinematic constraints restrict to be greater than the threshold energy squared 
for the process and the true invariant mass, M^^^^ = X1X2S = Q^/ {Z1Z2), to be greater 
than the visible invariant mass, M^. These result in the following constraints for Q^: 

Q' > Qthreshold , Q'>ZlZ2M\ (3.10) 

3.1 Single-invisible decays 

The benchmark scenario for a single invisible decay originating from the hard pro- 
cess is ti production in which one of the two tops decays into bqq' (hadronic) and the 
other into biu (semi-leptonic), the neutrino comprising the missing four- momentum. 
Excluding the proton remnants, we assume that all other particles within the pseudo- 
rapidity coverage are detected. We will refer to the neutrino as the invisible particle 
and the W as the intermediate particle in the tt case, but the treatment is readily ap- 
plicable to the gluino case where the invisible particle is the Xi and the intermediate 
particle is a squark (treated in section [3l^ ) . 

To calculate the functions and obtain Q^, we need to calculate the neutrino 
four-momentum in the hard process frame. This is done by choosing the neutrino 
four-momentum in the frame of its parent W and then applying two subsequent 
Lorentz boosts: one going from the W frame to the top frame, and one from the top 
frame to the hard process frame. The decay chain is shown in Fig. |3.1| . Each of these 
boosts involves two angular variables which originate from the 'decay' of the parent 
particle. Hence the four momentum of the neutrino may be written as 

gL = K {Q, 0, 0) Aa {0, 4>) PtiO, 0) , (3.11) 

where the A's are Lorentz boost matrices and where quantities with a hat refer to 
the hard process frame, quantities with a tilde refer to the top frame and quantities 
with a bar refer to the W frame. The angles 6 and represent the usual polar 



angles, defined with respect to the direction of the 'sister' particle (see Fig. |3.1|) . For 
example, in the case — )■ i'^i'i, where the was produced from the top decay 
along with a bottom quark, the angles {9, 0) are defined with respect to the direction 
of motion of the b in the frame. The two boost vectors have magnitudes given 
by = \pi\/Ei {i = t, W), the ratio of the parent 3-momentum magnitude and its 
energy. The boosts, as well as the magnitude of the invisible particle four-momentum, 
can obtained by considering kinematics in each frame as: 

pt{eA) = ^{l,r^, (3.12) 




Figure 2: The sequential two-body decay chain under consideration in the invisible particle 
treatment. The relevant production angles in the parent centre-of-mass frame are also 
shown in parentheses. 

where r = (sin 6 cos 0, sin 6 sin 0, cos 0) is the unit vector in spherical polar coor- 
dinates in the appropriate frame and m\Y, rrit are the W and top quark masses 
respectively. The four- vector f^!^^, and hence the functions f^^, are calculated by 
finv — Qinv/Q- Evidently, the functions are functions of Q^, giving an imphcit 
equation for Q^. To make this more explicit, we re- write eq. (p.6|): 



r)2 ^ ziZ2 [M^-ml^{Q^n)] 



and analogously for eq. ( |3.9| ), where Q represents the set of all angular variables. In 
the present case Q) = rriu ~ but for multiple invisible particles it will also 

be a function as indicated. 

Equation (|3.13|) needs to be solved numerically for each set (2:1,^2,^) in the 
region {Am'f^-, Z1Z2S), where S is the square of the proton centre-of-mass energy, 
along with the restriction that the visible invariant mass should be lower than the 
'true' invariant mass, M < Mtme- The numerical solution was found using the Van 
Wijngaarden-Dekker-Brent method [4,5], a bracketing method for finding roots of 
one- dimensional equations. Since Q is not uniquely determined for each M, dif- 
ferent values of the 'true' centre-of-mass energy Q contribute to the cross section. 
Note that not all possible configurations {zi, Z2, ^) are kinematically allowed to con- 
tribute to the cross section at M and hence some configurations do not yield roots 
of eq. ( |3.13|) . Once is obtained, the parton- level cross section for the hard pro- 



cess partons, o'a'b'iQ'^), is calculated. This result is then multiplied with the parton 
density functions for the incoming partons, fa,b{xi^2,Qc), and the kernels for evo- 
lution from incoming partons a and b to hard process partons a' and b' {Ka'a{zi) 
and Kb/b{z2)). We then integrate over all possible values of Zi and Z2, according to 
eq. ( p.7| ). Finally, to obtain the full resummed result we have to integrate over the 
distribution of the angular variables Q. Notice that the visible invariant mass distri- 
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bution becomes non-zero below the threshold for production, M < 2mt/g, owing to 
the loss of invisible particles. 

3.2 Double-invisible decays 

We now turn to the case where both particles produced in the hard process decay 
invisibly. For illustration we refer to sequential decays of the gluino: ^ — )■ gg — )■ XiQQ- 
Although this decay mode is generally not the dominant one, it is useful for illustra- 
tion of the procedure. We extend the treatment given in the semi-leptonic/hadronic 
top case by writing out functions related to the two invisible particle four-momenta 
in the decay chain (which we call x ^i^nd x'); 

= a;: (q, e, 0) A^ {e, 0) p'^ie, 0) , (3.i4) 

q^-, = A^ (q, 9', 0') A^ (O', 0') p'^,{e', 00 , (3.15) 

where the primed quantities now distinguish between the two invisibles. Since both 
of these four-vectors are defined in the hard subprocess frame, we have simply 

fL = ^{Qi + Qi) ■ (3.16) 

The rest of the treatment is identical to the one-invisible case: an implicit equation 
has to be solved to obtain for each [zi, Z2, ^) set and then an integral over Q is 
taken to obtain the resummed result. 

3.3 Angular distributions 

The distributions of the angular variables Q = {9, (j), 9, (p, 9,(1)), appearing in the 
treatment of invisibles given in the previous sections, are pro cess- dependent. They 
represent the angles at which the daughter particle is emitted in the frame of the 
parent particle. We investigated the angular distributions using a Monte Carlo event 
generator (Herwig++ 2.4.0 [3]) and subsequently used the results in calculating the 
functions. The results for SPSla gluino pair-production are shown in Fig. ^ 
where the uniform distributions are shown for comparison (red horizontal line). Fig- 
ure I shows the distributions as obtained for top pair-production. The neutrino angle 
in the W frame is also compared to the analytic calculation. As expected, all the 
angles, in both cases, were found to be uniform (not shown). The form of all the 
distributions can be justified using general spin considerations: 

9i: The angular distribution of the angle 9i at which the fermions are produced in 
the hard process frame is expected to have the form ~ 1 + /3 cos^ 9i, where /3 is 
a pro cess- dependent constant. 
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Oi. The angle 6*^, is defined between the direction of the daughter boson {W oi q) 
with respect to the direction of polarization of the parent [t or g) polarization. 
The angular distribution for a spin- up fermion parent is then given by [6]: 



dA^t 



-^t d cos 6i 



Pai cos Oi) 



(3.17) 



where ai is a constant and P is the modulus of the polarization of the parent. 
Since the production processes for both ti and gg are parity conserving, there 
is also an equal spin-down (A^^) contribution to the total distribution with the 
sign of ttj reversed. This results in a uniform distribution for cos6'j. 
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Figure 3: Monte Carlo results for the gluino pair-production decay chain angles. From 
left to right: the production angle of the gluino in the hard process frame, the angle of the 
outgoing squark in the gluino frame and the angle of the outgoing neutralino in the squark 
frame. The uniform distributions are shown for comparison. 



Oi. In gluino pair- production, the decay products of the squark, g, which is a scalar, 
are uniformly distributed in cos 9. In top pair-production, on the other hand, 
the decay W — )■ Ivi is parity- violating and the distribution of cos 9 is forward- 
backward asymmetric in the W frame [7]. The angle 9 (sometimes called \E', 
see e.g. [8]) is used experimentally to infer helicity information on the W . The 
distribution may be written as 

(3.18) 

where F^, Fr and Fq are the probabilities for left-handed, right-handed and 
longitudinal helicities of the W in top quark decay respectively. The SM pre- 
dictions, (Fj^, Fr, Fq) = (0.304, 0.001, 0.695), yield the curve shown on the right 
in Fig. |. 



1 diV _ 3 
iVdcos^ ~ 2 
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Figure 4: Monte Carlo results for the top pair-production decay chain angles. Prom left to 
right: the production angle of the top in the hard process frame, the angle of the outgoing 
W boson in the top frame and the angle of the outgoing neutrino in the W frame. The 
uniform distributions are shown for comparison. The neutrino angle in the W frame is also 
compared to the analytic calculation. 



The spins of the two produced fermions (tops or gluinos) are correlated and this 
may cause a degree of correlation between the distributions of particles in the decay 
chains. We investigated whether these correlations play an important role in the 
calculation of the invisible particle effects on the visible mass. By comparing the 
invariant mass distributions with and without the spin correlations in the Monte 
Carlo we concluded that the effect is small in both top and gluino pair-production 
and can be safely neglected. 



4. Results 

We present the resummed distributions obtained for tt and gg production according 
to eq. (|3.7] ). All results are for the LHC at design energy, i.e. pp collisions at a/s = 14 
TeV. We have integrated over the visible system rapidity, Y, in the range \Y\ < 5. 
We first compare our results to those obtained using the Herwig++ event generator 
at parton level (i.e. no hadronization or underlying event) and excluding the proton 
remnants.^ In sections and 14. 41 we examine the effects of hadronization and the 



underlying event. Parton- level top and gluino pair-production cross section formulae 
are given in appendix 0. The PDF set used both in the calculation and Herwig++ is 
the MRST LO** (MRSTMCal) set [9,10]. 

4.1 Top quark pair production 

We present resummed results in comparison to Monte Carlo for Standard Model ti 
production, where we include particles with maximum pseudorapidity rj^^x = 5. In 



^We verified using the event generator that the contribution of the proton remnants to the total 
invariant mass in the considered rapidity range is negligible. 
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Figure 5: The visible mass distributions for a pseudorapidity cut 7?max = 5. Left: com- 
paring hadronic (no invisibles) and semi-leptonic (one invisible) decays. Right: comparing 
hadronic (no invisibles) and fully leptonic (two invisibles) decays. The leading order tt 
invariant mass distribution is shown (red dot dashes) for comparison. 



Fig. ^ we show separate results for combinations of hadronic and semi-leptonic decays 
of the top, leading to zero, one or two invisible neutrinos from the hard process. The 
effect of the invisibles in both the fully semi-leptonic case and the hadronic/semi- 
leptonic case are small compared to the effects of hadronization, to be discussed in 
section O. The differences between the Monte Carlo and resummed curves in Fig. ^ 
may be attributed to sensitivity to the behaviour of the PDFs and parton showering 
at low scales, and the precise definition of Qc in terms of r^max, since Qc can be as 
low as 21711 X ~ 2 GeV in the case of ti production. 



4.2 Gluino pair production 

We focus on the SPSla point [11], which has gluino and lightest neutralino masses 
rrig = 604.5 GeV and m^o = 97.0 GeV respectively (and see table ^ for the squark 
masses). For simplicity we set the squark mass in the invisible particle treatment to 
550 GeV. We also present results for a modified SPSla point, with rrig = 800 GeV. 
In this process only the two-invisibles case is realistic, but for comparison we also 
show results for no invisibles, i.e. imagining that the two lightest neutralinos are also 
detected. When r/max = 5, 3, there is fairly good agreement between the Monte Carlo 
and resummation predictions in both the two-invisibles and no-invisibles cases, and 
for both gluino masses, as shown in Figs. ^ and |^ where one should compare the 
dashed histograms (Monte Carlo) to the solid curves of the same colour (resumma- 
tion). 

The shift in the peak of the visible mass distribution in going from no to two 
invisibles is much larger than that in top pair production, amounting to 600-700 
GeV, roughly independent of r^max and the gluino mass. This results mainly from 
the higher masses of the intermediate particles in the decays {rrig ~ 550 GeV vs. 
m\Y = 80 GeV), which implies a higher energy release, rather than the masses of the 
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Particle 


Mass (GeV) 


Particle 


Mass (GeV) 


9 


604.5 


Sl 


570.7 


X? 


97.0 


Sr 


547.9 


Ul 


562.3 


hi 


515.3 


ur 


548.2 


b2 


547.7 


dL 


570.7 


h 


400.7 


dR 


547.9 


t2 


586.3 



Table 1: The relevant particle masses in the supersymmetric model used in the invisible 
study, SPSla. The modified SPSla point differs in that it has rrig = 800 GeV. 



invisible particles themselves {fn^o = 97 GeV vs. = 0). 



One of the assumptions of the resummation is that all the visible hard process 
decay products are detected, which is not true when the maximum pseudorapidity 
'7max is restricted to lower values. When ?7max ~ 2 in the Monte Carlo analysis, 
a significant number of hard process particles begin to be excluded and hence the 
curves shift to lower values than the resummed predictions. Figure ^ shows the 
rapidity distribution of the decay products of the gluino at parton level for nig = 
604.5 GeV. For the case shown, cuts of rjma.x = 5, 3, 2 and 1.4 correspond to exclusion 
of, respectively, ~0.002%, 1.1%, 7.5% and 20.0% of the gluino decay products from 
the detector. The effect of this appears in Figs. ^ and |10], where the Monte Carlo 



distributions are narrower and peak at lower masses than the resummed predictions. 
The variation between the resummed ?7max = 2 and 1.4 curves is smaller than that 
between ?7jnax = 5 and 3, since they correspond to smaller differences in Qc- 

The heavy and light gluino scenarios exhibit similar behaviour when varying the 
pseudorapidity coverage and the number of invisibles, showing the lack of dependence 
of the resummation on the mass of the pair-produced particle. The sensitivity to low- 
scale PDF behaviour and showering is reduced compared to the tt case since we are 
considering higher centre-of-mass energies, with the lowest possible Qc now being of 
the the order 2m^ x ~ 8 GeV. The position of the curves is also sensitive to the 
precise definition of Qc in terms of ?7max- 

Table ^ shows a summary of the peak positions for all cases and different pseudo- 
rapidity cuts. For the higher values of ?7max; the agreement between the Monte Carlo 
and resummation is satisfactory. There is a large difference in the peak positions for 
no invisibles and r/max = 5, but this is mainly due to the broad shape of the peak 
in this case, while the overall distributions agree better. For ?7max < 2 there is a 
growing discrepancy, especially for the realistic case of two invisibles, due to the loss 
of particles coming from the hard process. 
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nig (GeV). 




MC (GeV) (0 inv./2 inv.) 


Resum. (GeV) (0 inv./2 inv.) 


604.5 


5 


2280/1560 


1785/1620 


604.5 


3 


1680/1080 


1593/1204 


604.5 


2 


1440/840 


1497/1204 


604.5 


1.4 


1380/660 


1497/1204 


800.0 


5 


2820/2100 


2569/1870 


800.0 


3 


2220/1620 


2128/1684 


800.0 


2 


1920/1380 


1865/1683 


800.0 


1.4 


1740/1140 


1865/1683 



Table 2: Summary of the positions of the peaks of the gluino pair-production visible mass 
distributions as given by the Monte Carlo and the resummation, for different values of the 
maximum pseudorapidity and for no and two invisibles. 



4.3 Hadronization effects 

We have assumed that ISR partons emitted at pseudorapidities above ?7max do not 
contribute to the visible invariant mass. This would be true if the hadronization 
process were perfectly local in angle. However, as a result of hadronization high- 
rapidity ISR partons can produce lower-rapidity hadrons and thus 'contaminate' the 
detector and shift the visible mass to higher values. 

The hadronization model employed in the Herwig++ Monte Carlo is a refinement 
of the cluster model described in ref. [12]. The model involves clustering of partons 
into colour-singlet objects that decay into hadrons, resulting in a smearing of the 
pseudorapidity distribution which causes the increase in the visible mass described 



above. The effect is shown in Fig. ITT] for gluino and top pair-production (excluding 
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M/GeV M/GeV 

Figure 6: The SPSla gluino pair-production visible mass distributions for pseudorapidity 
cuts r/max = 5 (left) and r/max = 3 (right). The leading order distribution is shown (red dot 
dashes) for comparison. 
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Figure 7: The modified SPSla gluino pair-production (with rrig = 800 GeV) results for 
pseudorapity cuts ?7max = 5 (left) and r/max = 3 (right) . The leading order distribution is 
shown (red) for comparison. 

the invisible particles from the hard process). The effect was found to be larger for 
tt production where the mass distribution is shifted significantly, whereas in gluino 
pair production the shift is negligible. 

4.4 Underlying event 

The underlying event, which is thought to arise from multiple soft interactions of 
spectator partons, is a further source of non-perturbative contributions to the visible 
mass. If represents the "hard" visible 4-momentum studied in earlier sections 
and represents that due to the underlying event, the total visible mass is given 
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Figure 8: The SPSla gluino pair-production pseudorapidity distribution for rrig = 
604.5 GeV. 
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Figure 9: The SPSla gluino pair-production results for pseudorapidity cuts r/max = 2 
(left) and r?max = 1-4 (right). The leading order distribution is shown (red dot dashes) for 
comparison. 



by 

= {Ph + Pu? = Ml + Ml + 2{EhEu - P.hPzu) 

= Ml + Ml + 2Mu^Ml+$^cos\i{YH - Yu) . (4.1) 

where we neglect transverse momentum associated with the underlying event. Thus, 
even if the visible invariant mass due to the underlying event is small, its effect on 
the overall visible mass may be enhanced through the last term on the right-hand 
side. 

The underlying event is simulated in Herwig++ by a multiple parton interaction 
model along the lines of ref. [13]. In this model, for the rapidity ranges considered 
here, the underlying event is approximately process-independent and exhibits little 




Figure 10: The modified SPSla gluino pair-production (with rrig = 800 GeV) results for 
pseudorapity cuts r?max = 2 (left) and r/max = 1-4 (right). The leading order distribution is 
shown (red) for comparison. 
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Figure 11: The tt fully semi-leptonic (left) and SPSla gluino pair-production (right, with 
nig = 604.5 GeV) visible mass distributions for a pseudorapity cut r/max = 5 with and 
without hadronization (black and red respectively). 



correlation with the rest of the event. Therefore, to a good approximation, the 
distributions of the variables related to the underlying event, Yu and Mu, can be 
determined once and for all at each collider energy. The process-dependence comes 
primarily through the dependence on Yh and Mh, which can be calculated using the 
resummation formula given in eq. ( p.7| ). The overall visible mass distribution can 



then be obtained by convolution using eq. ( 4.1 



The effects of including the underlying event in the visible mass distribution are 



shown in Figs. and |T3|for ti and gluino pair production, respectively. The multiple 
parton interactions push the peak value to substantially higher masses. The shift 
amounts to about 250 GeV at rj^ax = 3 and 1.2 TeV at rj^ax = 5, and is roughly 
process-independent. However, since the underlying event is approximately uncorre- 
lated with the hard process, the visible mass distributions can be reconstructed well 
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Figure 12: The tt fully hadronic visible mass distributions for pseudorapity cuts rjmax = 5 
(left) and r]max = 3 (right), with and without multiple parton interactions (black and red 
respectively) and the reconstructed curves (blue dot dashes). The r/max = 5 curve was 
reconstructed using the resummed results for the visible mass and rapidity, whereas the 
^max = 3 curve was reconstructed using the Monte Carlo visible mass and rapidity. 



- 16 - 



no mpi 
with mpi 

with mpi, reconstructed 





with mpi 

with mpi, reconstructed 



Figure 13: The SPSla gluino pair-production (with nig = 604.5 GeV) visible mass dis- 
tributions for pseudorapity cuts rjmax = 5 (left) and rjmax = 3 (right), with and without 
multiple parton interactions (black and red respectively) and the reconstructed curves from 
the Monte Carlo visible masses and rapidities (blue dot dashes). 



by the convolution procedure outlined above, as shown by the blue dot-dashed curves 



in Figs. and |T3|. These features of the underlying event will need to be validated 
by LHC data on a variety of processes. Accurate modelling of the underlying event 
is important for practically all aspects of hadron collider physics. 



5. Conclusions 



In this paper we have presented detailed predictions on the total invariant mass M 
of the final-state particles registered in a detector, as a function of its pseudorapidity 
coverage r/max- This quantity provides the dominant contribution to many global 
inclusive observables such as the new variable sj^^^ in ( |1.2| ), which can provide infor- 
mation on the energy scales of hard processes. We have extended the resummation 
method presented in [1] to include the effects of invisible particle emission from the 
hard process. We have considered the case of one or two invisible particles and pre- 
sented results for Standard Model top quark pair production and SPSla gluino pair 
production, obtained using a numerical Mellin moment inversion method. 

In the case of ti production the invisible particles are neutrinos from W-boson 
decays and their effect on the visible invariant mass distribution is small, even when 
both decays are leptonic. This is mainly a consequence of the small W-boson mass 
compared to the overall invariant mass, rather than the negligible neutrino mass. 
For gluino pair production the invisibles are a pair of massive LSPs from squark 
decays. The LSP mass is again small compared to the overall invariant mass, but 
the squark masses are not, leading to a substantial downward shift in the visible mass 
distribution, of the order of the squark mass. In both cases the resummed predictions 
are in fair agreement with Monte Carlo estimates of the position of the peak in 
the distribution, provided the pseudorapidity range covered by the detector is large 
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enough (r^max ^ 3). For ?7niax ~ 3, the difference between the Monte Carlo prediction 
and resummed predictions is of the order of 100 GeV for both the heavy and hght 
gluino SPSla points. The agreement becomes worse when the pseudorapidity range 
is restricted, due to particle loss from the hard process. Table |^ shows the positions 
of the peaks of the distributions for the Monte Carlo results from Herwig++ and the 
resummation. 

These comparisons were made with Monte Carlo visible mass distributions at 
parton level. We found that non-perturbative effects, especially the underlying event, 
tend to shift the invariant mass distributions to significantly higher values than 
expected from a purely perturbative calculation. According to the underlying event 
model used in Herwig++, the shift amounts to about 250 GeV at rjmax = 3 and 
1.2 TeV at rj^ax = 5. This effect is also expected in other observables sensitive 

1/2 

to longitudinal momentum components, such as s^^jj^. However, in this model the 
underlying event is only weakly correlated with the rest of the event and hence its 
effects can be determined once and for all at each collider energy. The modelling of 
the underlying event is an important feature of the Monte Carlos that needs to be 
validated by comparison with experiment. Once this has been done, a wide range 
of global inclusive observables, including the visible invariant mass, will be reliably 
predicted and useful for establishing the scales of contributing hard subprocesses. 
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A. Calculation of the evolution kernels 

The evolution kernels, Ka'a{z), were calculated using the Mellin inverse transform 
method. Recall that the Mellin transform is defined by the expression 

f dzz'^-^Ky.iz) (A.l) 
Jo 

with inverse 

Ka'a{z) = ^ / <lNz-''Ki^ , (A.2) 

where the contour C runs parallel to the imaginary axis and is to the right of all 
singularities of the integrand. 
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Taking into account the running of the strong couphng, the Melhn transform of 



the solution of the evolution equation (|2.4| ) is 



K 



ba 

N 



as(Q) 



(A.3) 



ba 



with p = 6/(1 ICa — 2n/) and 



TT 



{^N)ba — iXN)ba 

as 



Jo 



(A.4) 



When computing this expression we need to take into account the fact that Ajv is a 
matrix. This can be done most easily by writing it as 



asiQc) 
asiQ) 



-1 pdiag(Ajv,i) 



o 



(A.5) 



where O is the matrix of eigenvectors of A^v and diag(A7v,'^) is the diagonal matrix 
of its eigenvalues. This is equivalent to using, implicitly, the singlet and non-singlet 
basis for the PDF evolution [7]. 
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Figure 14: Integration contours, Cq and Ci, for the inverse Mellin transform as given by 
the Bromwich integral, eq. ( [A.2| ). 

To evaluate the inverse transform (|A.2| ), numerical convergence of an otherwise 
infinitely oscillatory expression is achieved by choosing a contour that introduces a 
damping factor. A method commonly used in PDF evolution (see, for example, [14]) 
is to rotate the upper and lower portions of the vertical contour so that they slope 
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back into the left half-plane, as shown by Co in Fig. 14. This introduces exponential 
damping along the contour, which is chosen so as to enclose the same singularities as 
the 'Bromwich' contour C, and therefore converges to the correct result by Cauchy's 
theorem. However, in the case of the evolution kernels Ka>a{z) the linear contour Cq 
does not provide sufficient accuracy to reproduce the function from its transform. 
This is due to its inability to invert the constant function /^v = c to the correct 
analytic result, a delta function. This implies that the inversion does not reproduce 
the necessary initial condition of the evolution, Ka'a{z,Q = Qc) cc 5{1 — z). A 
numerically more accurate contour is available in the literature, used in the so-called 
'Fixed- Talbot algorithm'. This contour Ci has the form Re(iV) = ImA^cot(ImA^/r) 
where r is a parameter which we set to r = 0.4m/ log(l/z) during the computation, 
m being the required precision in number of decimal digits, a value derived from 
numerical experiments. The contour is related to the steepest descent path for a 
certain class of functions. For further details on its origin and accuracy see [15]. 



B. Pair-production cross sections 



The leading-order parton-level cross section for QCD pair-production of particles of 
mass nip may be written in terms of scaling functions fij as 



CTijiQ ) = — /) 



(B.l) 



For heavy quark pair-production, the functions for gluon-gluon and quark-antiquark 
initial states are given by [7]: 



7r/3p . , 



27 
' 192 



■i(p2 + 16p + 16) log 



1 + /3 
1-/3 



-28-31p 



(B.2) 
(B.3) 



where p = im'^/Q'^ and (3 = y/1 — p. For the case of gluino pair-production, the 
equivalent functions fij are given by [16]: 
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fqq 



+ 



9 9m? 
- H 

4 g2 
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9m^ 



log 

2 



l + (3 
1-/3 



3/3 
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4g2 



6m 



27 



+ 



32m^ 



64m| 8m^ 



16m? ; 

"gg^ 3Q2 ' 27(mi+m?g2)_ 



(B.4) 
(B.5) 



16m~m_ 



_27g2 3g4 27g2(g2-2m2 



log 



1 + /3 - 2m^/g^ 



where now /3 



1 — 4m?/g2 and represents the mass-squared difference be- 



tween the gluino and the t-channel squark, m_ 



2 2 

m| — m|. 



- 20 - 



References 



A. Papaefstathiou and B. Webber, "Effects of QCD radiation on inclusive variables 
for determining the scale of new physics at hadron colliders," JHEP 06 (2009) 069, 



0903.2013. 



P. Konar, K. Kong, and K. T. Matchev, ^^yf^rain ■ A Global inclusive variable for 
determining the mass scale of new physics in events with missing energy at hadron 
colliders," JHEP 03 (2009) 085, |0812 . 1042| . 

M. Bahr et al, "Herwig++ Physics and Manual," Eur. Phys. J. C58 (2008) 



639-707, 0803.0883 



R. Brent, Algorithms for Minimization without Derivatives. NJ: Prentice-Hall, 1973. 

M. Galassi, K. Davies, J. Theiler, J. Cough, G. Jungman, P. Aiken, M. Booth, and 
F. Rossi, GNU Scientific Library Reference Manual (3rd Ed.). Network Theory Ltd, 
2009. 

F. Hubaut, E. Monnier, P. Pralavorio, K. Smolek, and V. Simak, "ATLAS sensitivity 
to top quark and W boson polarization in ti events," Eur. Phys. J. C44S2 (2005) 



13-33, |hep-ex/0508061 



R. Ellis, W. Stirling, and B. Webber, QCD and Collider Physics. Cambridge 
University Press, 1996. 

The ATLAS Collaboration, C. Aad et al, "Expected Performance of the ATLAS 
Experiment - Detector, Trigger and Physics," [0901 . 0512 . 



A. Sherstnev and R. S. Thorne, "Parton Distributions for LO Generators," Eur. 



Phys. J. C55 (2008) 553-575, |0711.2473 



A. Sherstnev and R. S. Thorne, "Different PDF approximations useful for LO Monte 
Carlo generators," |0807.213^ . 



B. C. Allanach et al, "The Snowmass points and slopes: Benchmarks for SUSY 
searches," Eur. Phys. J. C25 (2002) 113-123, |hep-ph/020"223^ . 



B. R. Webber, "A QCD Model for Jet Fragmentation Including Soft Gluon 
Interference," Nucl. Phys. B238 (1984) 492. 

J. M. Butterworth, J. R. Forshaw, and M. H. Seymour, "Multiparton interactions in 



photoproduction at HERA," Z. Phys. C72 (1996) 637-646, |hep-ph/9601371 . 



A. Vogt, "Efficient evolution of unpolarized and polarized parton distributions with 



QCD-PEGASUS," Comput. Phys. Commun. 170 (2005) 65-92, tiep-ph/0408244. 



J. Abate and P. Valko, "Multi-precision Laplace transform inversion," International 
Journal for Numerical Methods in Engineering 60 (2004) 979-993. 



- 21 ~ 



[16] W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas, "Gluino pair production at 



the Tevatron," Z. Phys. C69 (1995) 163-166, |hep-ph/9505416 



- 22 - 



